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ABSTRACT 

Empirical 'backward' galaxy evolution models for IR-bright galaxies are constrained using 
multiband IR surveys. A new Monte-Carlo algorithm is developed for this task. It exploits a 
large library of realistic Spectral Energy Distributions (SEDs) of 837 local IR galaxies (IRAS 
25iim selected) from the UV (1000A) to the radio (20cm), including ISO-measured 3-13/im 
unidentified broad features (UIBs). The basic assumption is that the local correlation between 
SEDs and Mid-Infrared (MIR) luminosities can be applied to earlier epochs of the Universe, 
an assumption which will be strongly tested by SIRTF. By attaching an SED appropriately 
drawn from the SED library to every source predicted by a given model, the algorithm 
enables simultaneous comparisons with multiple surveys in a wide range of wavebands. Three 
populations of IR sources are considered in the evolution models. These include (1) starburst 
galaxies, (2) normal late-type galaxies, and (3) galaxies with AGN. Constrained by data from 
the literature, our best-fit model ('Peak Model') predicts that since z=1.5 the population of 
starburst galaxies undergoes a very strong luminosity evolution (L = Lq x (1 + z) 42 ) and also 
strong density evolution (p = po x (1 + z) 2 ), the normal late-type galaxy population undergoes 
a passive luminosity evolution (L = L x (1 + z) 1 - 5 ), and the galaxies with an AGN undergo a 
pure luminosity evolution similar to that of optical QSOs (L = L x (1 + z) 35 ). Prior at z>1.5 
all evolution rates drop as (1 + z)~ 3 . The luminosity evolution results in evolution of SEDs of 
IR-bright sources because of the luminosity dependence of the SEDs. Predictions for number 
counts, confusion limits, redshift distributions, and color-color diagrams are made for multiband 
surveys using the upcoming SIRTF satellite. A A-Cosmology (fl\ = 0.7, il m = 0.3, Ho=75 km 
sec -1 Mpc -1 ) is assumed throughout the paper. 

Subject headings: galaxies: starburst - Seyfert - luminosity function; infrared: galaxies 



1. Introduction 

The first sign of cosmic evolution among infrared (IR) galaxies was detected by Hacking, Condon and 
Houck (1987) in the IRAS 60/ito deep survey (Hacking and Houck 1987). This was subsequently confirmed 
by later studies of IRAS galaxy populations (Franceschini et al. 1988; Lonsdale and Hacking 1989; Lonsdale 
et al. 1990; Rowan-Robinson ct al. 1990; Saunders et al. 1990; Yahil ct al. 1991; Grcgorich ct al. 1995; 
Pearson and Rowan- Robinson 1996; Bertin et al. 1997). Recently, deep MIR-FIR surveys have been carried 
out using the Infrared Space Observatory (ISO) (Kessler et al. 1996). These include ISOCAM surveys at 
15iim, 12/im and at 6.7/xrn (see Elbaz et al. 1998b for a summary of these observations), ISOPHOT surveys 
at 90/xm (Oliver et al. 2000, Efstathiou ct al. 2000a) and at 175^m (Kawara et al. 1998; Puget et al. 1999; 
Dole et al. 2001). The results from these surveys (Aussel et al. 1999; Puget et al. 1999; Dole et al. 2001; 
Clements et al. 1999; Elbaz et al. 1999; Serjeant et al. 2000; Xu 2000) indicate strong cosmic evolution in 
the population of infrared-emitting galaxies, confirming the earlier results based on smaller samples and 
less sophisticated analyses (e.g. Rowan- Robinson et al. 1997; Kawara et al. 1998). This is consistent with 
the results of SCUBA surveys (Hughes et al. 1998; Barger et al. 1998; Blain et al. 1999) and with the 
scenario hinted at by the newly discovered Cosmic IR Background (CIB) (Puget et al. 1996; Hauser et al. 
1998; Dwek et al. 1998; Fixsen et al. 1998), while challenging the results from UV/optical surveys in the 
sense that substantially more (i.e. a factor of 3 - 5) star formation in the earlier Universe is required to 
match the IR/submm counts and the CIB (see, e.g. Rowan-Robinson et al. 1997) compared to that derived 
from the UV/optical surveys (Madau ct al. 1998; Pozzetti et al. 1998, Stcidcl et al. 1999). The reason for 
this discrepancy is attributed to dust extinction which may hide much of the star formation in the early 
Universe from the UV/optical surveys (see Lonsdale 2000 for a review). 

Compared to the UV and optical surveys, the infrared surveys are superior in their insensitivity to 
dust extinction, but are inferior in angular resolution (a few arcsecond compared to the sub-arcsecond 
resolution of optical surveys). This not only limits the IR surveys by confusion, but also makes the study 
of the IR morphology of faint IR sources impossible. In order to reveal the true nature of faint IR sources, 
identifications in other bands, especially in optical and NIR bands where the sources can be resolved easily 
with current instruments, are usually needed. The multi-band studies (including ISOCAM 15/im and 
6.7/um surveys) of the Hubble Deep Field (HDF) by Rowan-Robinson et al. (1997) and Aussel et al. (1999), 
and of the Canada-France Redshift Survey (CFRS) field by Flores et al. (1999) suggest that, compared to 
their optical counterparts, the ISOCAM sources have significantly redder (I-K) colors (Flores et al. 1999), 
and are much more likely to be in the interacting/merging systems. On the other hand, such studies are 



necessarily confined to IR sources that are relatively bright in the optical (e.g. I < 22.5 mag, Flores et al. 
1999), while many IR bright galaxies are optically faint due to heavy dust extinction. 

IRAS studies showed that galaxies of different nature in the local universe have distinct IR spectral 
energy distributions (SEDs). Galaxies with bright active galactic nuclei (AGN) usually have significantly 
lower /60/i//25/i ratios (dc Grijp et al. 1985; Fang et al. 1998) than other galaxies. Interacting/starburst 
galaxies such as M82 have systematically high /60/i//ioo^ an d /25^//i2/i ratios than normal galaxies such as 
the Milky Way (e.g. Hclou 1986). In principle, these different characteristics in the IR SEDs for different 
populations of galaxies should facilitate a tool for identifications of IR galaxies when multi-band IR surveys 
are available, independent of the optical identifications. For ISO surveys this may not be very relevant 
because the ISOPHOT FIR surveys do not match the ISOCAM MIR surveys in depth due to severe 
degradation of the sensitivity of ISOPHOT detectors. However, when the Space Infrared Telescope Facility 
(SIRTF) is launched in mid 2002, simultaneous deep surveys in seven MIR-FIR bands (3.6/im, 4.5^m, 
5.8/im, 8.0/im, 24/im, 70/um, and 160/im) will be possible (Bicay et al. 1999). These will include the 
Guaranteed Time Observer programs with MIPS and IRAC (sirtf.caltech.edu/ROC/Titles_abstracts.html), 
the large-area Legacy survey SWIRE (Lonsdale et al. 2001; and the very deep Legacy survey GOODS 
(Dickinson et al. 2001). 

Significant K corrections will occur in the observed SEDs of faint sources in deep IR surveys. In 
particular, in the rest frame wavelength range between 3-20 /im there are several broad band features 
(see Puget and Leger 1989 for a review), often referred to the Unidentified Infrared Bands (UIBs), which 
are ubiquitously present in the MIR spectra of local galaxies with equivalent widths up to several microns 
(Hclou et al. 2000), with the exception of type 1 Seyferts (Clavcl et al. 2000). If these features are also 
present in the SEDs of high redshift galaxies, substantial K-corrections will result in when any of the 
features redshifts in or out of the band pass of an IR filter. These effects may indeed be beneficial rather 
than annoying, for they may facilitate IR photometric redshift techniques. 

Xu et al. (1998, hereafter Paper I) studied the effect of K-corrections due to UIBs on number counts 
of MIR surveys. In that work, a three-component model, with empirically determined MIR SED templates 
of (1) a cirrus/PDR component (2) a starburst component and (3) an AGN component, is developed for 
infrared (3-120 /um) SEDs of galaxies. The model is then applied to a complete IRAS 25 /im selected 
sample of 1406 local galaxies (z < 0.1; Shupe et al. 1998). Results based on these 1406 spectra show 
that the MIR emission features cause significant effects on the redshift dependence of the K-corrections, 



which in turn affect deep counts and rcdshift distributions in MIR surveys. In a subsequent paper, Xu 
(2000, hereafter Paper II) found that indeed the sharp peak at about 0.4mJy in the Euclidean normalized 
differential counts at 15/im (Elbaz et al. 1999) can be explained by the effects of UIBs, together with 
an evolution rate significantly stronger than derived in previous IRAS studies, eliminating the need for a 
hypothetical 'new population' (Elbaz 1998a). 

In this paper, we expand the models in Paper I in several aspects. 

(1) First of all, the analytical algorithm of the number count model which includes a proper treatment of 

the K-correction (Eqs (23), (24) and (25) in Paper I) is replaced by a Monte-Carlo algorithm, in which 
every source (galaxy) in a volume of given rcdshift and in a given luminosity bin is assigned an SED 
appropriately selected from the SED library (837 SEDs). The source's flux densities in different bands 
are then calculated by convolving the redshifted SED with the band passes of filters. In this way, we 
effectively simulate a virtual sky for a given evolution model. This not only enables the simultaneous 
comparison of counts in different bands, but also preserves the correlations between flux densities of 
different bands. This latter feature allowing us to predict color-color diagrams of different populations 
as a function of redshift, facilitating the exploration of photometric redshift indicators. 

(2) In the 'backward evolution' model, instead of treating all IR sources as a single population, in this 

work they are separated into three populations, in a similar spirit as in the model of Franceschini ct 
al. (1988; see also Roche and Eales 1999): (1) normal late-type galaxies, (2) interacting/starburst 
galaxies, and (3) galaxies with AGNs. These different populations are assumed to have different 
cosmic evolution rates. 

(3) The wavelength coverage of our SED library is expanded from 3 — 120/ito to lOOOA — 20cm. This is 

done by collecting from the literature the optical/NIR (B, J, H, Ks bands) magnitudes and the radio 
continuum (20cm) flux densities for galaxies in our SED sample, and by extrapolating from IRAS 
60/ito and lOOfim bands to submm bands using empirically determined correlations. 

The goal of this paper is to provide a set of comprehensive 'backward evolution' models (a category 
of galaxy evolution models in which number densities and other properties, e.g. luminosities in different 
bands, of local galaxies are evolved 'backward' in time from the present — i.e. with increasing redshift - 
according to some parametric prescriptions, see Lonsdale 2000 for a review) for future multi-band surveys, 
in particular those to be conducted with SIRTF. The model parameters will be constrained by observations 



available in the literature. This not only includes the ISO deep surveys, but also the optical/NIR surveys, 
SCUBA surveys and radio deep surveys, because our SEDs now cover all of these wavebands. Constraints 
derived from the Cosmic IR Background (CIB) will also be incorporated. The strength of models being 
constrained by such wide range of data has already been demonstrated in several previous papers (e.g. Blain 
et al. 1999; Trentham, Blain and Goldader 1998; Adelberger and Steidel 2000; Rowan-Robinson 2001). 

The model focuses on IR-bright galaxies, and therefore is not expected to match observed number 
counts in any band for which there is a substantial contribution from IR-quiet populations, eg. the K band, 
because E/S0 populations are missing from our model, and the bright radio counts, which are dominated 
by radio galaxies. 

Throughout the paper, the cosmology model specified by the following parameters is adopted: Ho=75 
km sec" 1 Mpc" 1 , Cl m = 0.3, fl A = 0.7. 

2. SED Library 

As in Paper I, our SED sample is drawn from the IRAS 25 /im selected sample (complete down to 
/25/i = 0.25 Jy) by Shupe et al. (1998), and contains 1455 galaxies, 1406 of them with redshifts < 0.1. As 
pointed out by Spinoglio et al. (1995), the MIR luminosities correlate well with the bolometric luminosities. 
Therefore MIR selected samples, such as ours, have fair representations of different populations of IR 
sources. However E/SOs, which are ~ 20% of optically selected galaxy samples but mostly undetected by 
IRAS, are not included in our SED sample. 

In Paper I, a three-component (cirrus/PDR, starburst, and AGN) MIR-SED model is applied to these 
galaxies, predicting an SED from 3-120 /im for each of them. In order to expand the SEDs to the optical and 
NIR bands, B (4400A), J (1.0 /im), H (1.6 /im), and Ks (2.2 /im), we searched the literature. B magnitudes 
of 1339 galaxies were found in the NASA/IPAC Extragalactic Database (NED). The NIR magnitudes are 
taken mainly from the 2MASS Second Incremental Data Release via the IRSA facilityn where J, H, and Ks 
magnitudes of 790 galaxies were found. In addition, NIR magnitudes of 413 galaxies in our sample are given 



1 The NASA/IPAC InfraRed Science Archive (IRSA) is a NASA project focused on providing software 
and Internet services to facilitate astronomical discoveries, support the production of new astronomical data 
products, and to plan future observations utilizing the data archives from infrared astrophysics missions 
supported at Infrared Processing and Analysis Center (IPAC). 



in Spinoglio ct al. (1995), 244 of which are overlapped with the 2MASS matchs. Whenever NIR magnitudes 
are available from both 2MASS and Spinoglio ct al. (1995), 2MASS data take precedence. Altogether, J H 
and Ks magnitudes are found for 959 galaxies in our sample. The radio continuum flux densities at 20 cm, 
Siaghz, were searched for in both the NVSS (Condon et al. 1998) and FIRST (Becker et al. 1995) catalogs 
of NRAO. Among 1406 galaxies in our sample 1170 are found in one of these two surveys (the rest are in 
the sky area not visible by VLA). It is found that 854 galaxies in our sample have B, J, H, Ks magnitudes 
and radio continuum flux at 1.4GHz. After excluding 17 galaxies which were undetected by IRAS both in 
the 60 /xm and 100 /im bands (whose IR SEDs are highly uncertain), we select a final SED sample of 837 
galaxies. Note that the 569 galaxies in the original sample (1406 galaxies) that do not make it into the final 
SED sample are mostly galaxies without NIR magnitudes. This is mainly due to the fact that data for NIR 
sources in a large fraction of the sky has not been releaseded by the 2MASS survey (the major source of 
the NIR data) yet. Since both the 2MASS survey and the VLA surveys are much deeper than the IRAS 
survey, the NIR magnitudes or the 20 cm flux density is missing for a source in the 25/zm selected sample 
only when the sky region is missing in the corresponding database. Therefore no bias is introduced into the 
final SED sample when sources without NIR or radio fluxes are excluded. 

For each galaxy in the SED sample, the broad-band UV-optical-NIR (1000A - 4 /um) SED is estimated 
by a spline fit of the fluxes in B, J, H, Ks bands, altogether with the predicted 4 fim flux density from the 
MIR SED model (Paper I). It should be noted that at the UV wavelengths (1000 — 4000A) the predicted 
fluxes are extrapolations from the available data, and caution should be applied when these predictions are 
used. This aspect of the model will be improved in the next paper, using the new UV data that are just 
now becoming available for ULIRGs and other IR-bright galaxies. 

The MIR (4 - 16/ito) SED is determined using the MIR SED model developed in Paper I, including a 
full treatment of the UIB features (absent in type 1 AGN). Then the SED in the wavelength range 16 — 
1200 [im is specified by a spline fit of IRAS data at 25, 60, 100 /im, together with the 16 /im flux density 
predicted by the MIR SED model, and the 170/xm, 240/ito, 450/ito, 850/ito and 1200/ito flux densities 
predicted by empirical correlations between the given submm band flux and the IRAS 60/ito and 100 /im 
fluxes which are derived from available submm data collected from the literature (Appendix). 

The radio continuum flux density at 20cm (Siaghz) is extrapolated to 6.2cm (5*4. sghz) and 2.8cm 
(Sio.2GHz) using the mean spectral indices a20cm/e.2cm = 0-79 and a6.2cm/2.8cm = 0.70, found for 
Shapley-Ames galaxies (Niklas et al. 1997). These radio flux densities are then linked to the end of the 



IR/submm SED at 1200^m by spline fit. 

In Fig.l we show the SEDs which are binned according to population (see the next section) and the 
25/i?n luminosity (Table 1). For each bin, the mean SED and its la dispersion are also plotted in the 
corresponding panel in Fig.l, and the values are listed in Table 2. 

As examples, in Fig. 2 observational data of twelve well known galaxies are compared with model SEDs. 
The data are collected from the literature, and the sources are 

• broad band optical magnitudes: NED; 

• NIR magnitudes: 2MASS and Spinoglio et al. (1995); 

• FIR/submm flux densities: IRAS, Bcnford (1999), Dunne et al. (2000), Rigopoulou et al. (1996), 
Lisenfeld et al. (2000), Carico et al. (1992), Chini et al. (1986), Andreani and Franccschini (1996), 
Roche and Chandler (1993); 

• radio continuum flux densities: Niklas et al. (1995), Condon et al. (1990). 

The agreements between the data and the model SEDs are remarkably good in general. 

3. 'Backward Evolution' Model for Multi-band Surveys 

3.1. Three Populations of IR Emitting Galaxies and Their LLFs 

In Paper I and Paper II, it is assumed that all IR sources evolve as a single population. Here we improve 
on that formulation by adopting a model in which IR sources can be separated into three populations, in 
a similar spirit as in the model of Franccschini ct al. (1988; see also Roche and Eales 1999): (1) normal 
late-type galaxies, (2) interacting/starburst galaxies, and (3) galaxies with AGNs. 

To enable these different galaxy populations to evolve at different rates, the model requires a local 
luminosity function (LLF) for each component. We began with the 25/im flux density-limited sample of 
1455 galaxies of Shupe et al. (1998). Although classifications of many of the galaxies as AGNs, normal 
galaxies, etc. are available in databases such as NED, to treat the sample in a more uniform way, we use 
IR-color-based criteria to divide our sample into different populations (see Paper I and Fang et al. 1998 for 
a discussion). We chose the following IRAS color boundaries: 



• "AGNs": f(j0fim/f25^m < 5 

• "Starbursts" : fea^m/ 32h\xm > 5 and /ioo^m/ fm^m < 2. 

• "Normals": feofim/ hs^m > 5 and fioo^m/ ' fm^m > 2. 

This division results in 356 galaxies classified as "AGNs" , 456 galaxies as "Normal" , and 643 galaxies as 
"Starburst" . We note that our color-selection method for AGN will not be appropriate for heavily obscured 
objects in which even the feOfim/ f25p,m color may be affected by reddening; instead any such IR-red AGN 
will be found in the 'starburst' sample. Also, for many AGNs as defined above, much of the IR radiation 
can be due to the emission of dust heated by stars in the host galaxy in addition to the dust emission 
associated with the AGN. For example, according to Eq(l) of Paper I, an IR source of feotim/ f25fim = 5 
(an 'AGN' by the above definition) has a half of its 25/xm emission from dust heated by stars. In reality, 
Seyfert galaxies such as NGC 4945 (Spoon et al. 2000) can have the IR emission predominantly powered 
by the nuclear starburst rather than by the AGN. This explains why the SEDs of some sources in the AGN 
subsample of our SED library show significant broad band MIR emission features (Fig. la) which should be 
absent in a typical type 1 AGN SED (e.g. Fig. 3 of Paper I). At the same time, as shown in Fig. 2, the UIB 
features are indeed absent in the model SEDs of AGNs such as Mrk 231 amd NGC 7479, consistent with 
the fact that the emission in these sources is dominated by the AGN. It should also be noted that in our 
SED model (Paper I), we do not distinguish type 1 and type 2 AGNs, which have significantly different 
MIR SEDs (Clavel et al. 2000). Many type 2 AGNs have strong UIB features even when their IR emission 
may be predominantly from dust associated with AGN, because an edge-on torus may heavily extinguish 
the MIR part of the AGN-associated emission and therefore the detected MIR emission is mostly from dust 
in the ISM of the host galaxy (Clavel et al. 2000). 

Luminosity functions were then computed for each of these populations according to the maximum- 
likelihood method described in Yahil et al. (1991) and used in Shupe et al. (1998). This method calculates 
the shape of a parametric luminosity function described by the parameters a, (3 and L* independent of 
density variations. 

With the shape parameters in hand, the normalization of each luminosity function must be estimated 
by other methods. Since the normalization of the total 25 micron luminosity function was estimated in 
Shupe et al. (1998), the normalizations of the three population LFs are chosen so that the the number of 
galaxies implied by the sum of the component LFs is about the same as the total LF. The difference between 
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the total LF and summed LF is also constrained to be less than a few percent at all luminosities. We have 
adjusted the normalizations of the component LFs to satisfy these criteria. These relative weightings of the 
component LFs can be adjusted by ten to twenty percent while still satisfying the criteria, but cannot be 
made vastly different from the nominal values. 

A plot of the component luminosity functions, the sum of the component LFs, and the total LF is 
shown in Fig. 3. The computed parameters for each population are given in Table 3. Note that the LFs 
are calculated using the whole sample (1455 galaxies) of Shupe et al. (1998). Identical results are obtained 
when the sample is confined to the 1406 galaxies with z < 0.1. 

3.2. Monte-Carlo Simulation of Multi-band Surveys 

In this subsection, we develop the algorithm to model coherently the number counts in different bands 
and the color-color diagrams for multi-band surveys. 

In a flat A-Universe (i.e. fi = fl m + fl\ = 1, Q\ =/= 0), which is adopted in this work, the co- moving 
volume is 

V=^D\ t (1) 

where A is the sky coverage in steradians, and Dm is the proper motion distance (Carroll et al. 1992)g: 

c f zl 
D m = tt [(1 + zf{l + Q m z) - z(2 + z)Q A ]-^ 2 d z (2) 

"0 Jo 

The predicted number of sources from a given population, in a given redshift interval [z — 0.5<5z, z + 0.5Sz] 
and in a given 25fmi luminosity interval [L — 0.5SL, L + 0.55L], is then 

6N i (L,z) = p' i (L,z)^-5L6z (3) 

where p\ is the luminosity function of population i (i = 1 - normal late type galaxies, i = 2 — starburst 
galaxies, and i = 3 - galaxies with AGNs): 

d(L,z) = G i (z)p i (j^y, (4) 



2 The Da, the angular diameter distance in Section 4.1 of Paper I, should have been called Dm, the 
proper motion distance, too. 
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where pi is the local 25 pxa luminosity function (luminosity function at z—0) of population i (Section 3.1), 
and Gi(z) and Fi(z) are the density evolution function and the luminosity evolution function of population 
i, respectively. 

For each of the SNi(L, z) sources predicted using Eq([f), an SED is randomly selected from a proper 
SED bin in the SED library (Section 2) and is assigned to the source. The SED library is binned according 
to (1) the population, and (2) the 25/im luminosity (Table 1). For an individual source the flux density in 
a given band can then be determined as follows: 

f band = tzW x L2 ld 2 lr I" 2 Si^-JR^^dX (5) 



AttD 2 l S{25pm) J Xl v l 

where 

D L = (1 + z) x D M (6) 

is the luminosity distance (Carroll ct al. 1992), Z^s^m — vL v (25pm) = XL \ (25/im) is the monochromatic 
luminosity at 25pm, Rband(X) is the bandpass of the given band, and S(X) is the flux density distribution 
of the SED in question. Due to the dependence of SED shape on luminosity for IR-bright galaxies (the 
/60/i//ioo/i color increases with Lju), this approach empirically results in color evolution accompanying 
luminosity evolution. 

When the evolution functions in Eq(4) are specified, we can predict counts in different IR bands, as 
well as contributions from IR galaxies to counts in other wavebands. As a test, local luminosity functions in 
the IRAS 60pm band (Fig.4), in the SCUBA 850 urn band (Fig.5), and in the IRAM 1250pm band (Fig.6) 
are calculated via model simulations (for sources of z < 0.1) specified by Gi{z) = 1 and Fi(z) = 1. Good 
agreement with the IRAS 60/ito luminosity function of Saunders et al. (1990), with the 850pm luminosity 
function of Dunne et al. (2000), and with the 1250pm luminosity function of Franceschini et al. (1998) is 
found. 

3.3. Evolution Models 

As in Paper I and Paper II, the following (power-law) function forms are adopted for the luminosity 
evolution functions Fi(z) and the density evolution functions Gi{z): 

F t (z) = (l + z) u * (z< Zl ) 

= (l + z) v * ( Zl <z<z ) (7) 
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Gi(z) = (l + z) Pi (z<zi) 

= (1 + z) qi (zi < z < 2d) (8) 

where zo is the redshift when the galaxy formation started, and z\ is the so called 'peak' redshift, where the 
evolution reaches a peak. Here we explore two kinds of models, the first characterized by a steady increase 
in evolution from z to Z\ at power law rates Vi and Qi followed by a strong decline to the present day 
with power law rates Ui and pi, and the second having a plateau between the formation and peak epochs, 
Zo > z > Z\. 

Throughout the paper we adopt zq — 7 and z\ = 1.5. Our results (predictions for number counts and 
the CIB) are not sensitive to zq so long as it is > 5. Optical surveys (e.g. Lilly et al. 1996; Madau et al. 
1996; Connolly et al. 1997) show that star formation rate in galaxies peaks between z = 1 and z = 2. The 
deep ISOCAM counts (e.g. Elbaz et al. 1999) are consistent with a rapidly increasing star formation rate 
going back at least to z > 1.5 (Paper II). The choice of 1.5 for the peak redshift is driven by the deep ISO 
15/iTO data (see Fig. 7). It is still controversial whether the star formation rate indeed decreases at z > Zi, 
beyond the peak, or flattens (e.g Steidel et al. 1999; Blain et al. 1999). 

Other parameters are Ui,Vi,Pi,qi (i=l, 2, 3). In order to reduce further the parameter space, we 
assume: 

(1) For normal late type galaxies (population 1): u\ = 1.5, p\ = 0. This corresponds to a pure passive 

luminosity evolution before the turn-over redshift. These galaxies are a major constituent of K-band 
cxtragalactic source counts, which may be consistent with passive evolution models (Gardener et al. 
1997). 

(2) For galaxies with AGNs (population 3): U3 = 3.5, pj, = 0. Here the assumption of the pure luminosity 

evolution is based on the studies of the evolution of QSOs in the literature (e.g. Boyle et al. 1988; Pei 
1995; La Franca et al. 2000). Using a maximum-likelihood technique Boyle et al. (1988) found that 
pure luminosity evolution models of form L(z) oc Lq x (1 + z) 1 (7 = 3.2 ± 0.1) adequately describe 
the evolution of bright (Mb < —23) low reshift (z < 2.2) QSOs. A similar result was found by Pei 
(1995), with 7 in the range of 3.2 — 3.9. Recently, La Franca et al. (2000) found that results from 
an ISOCAM 15/ito survey of Type 1 AGNs in the ELAIS fields are consistent with a pure luminosity 
evolution model with L(z) oc Lq x (1 + z) 3A . 

(3) The evolution indices ui and P2, which specify the luminosity and density evolution rate of starburst 
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galaxies before z = 1.5, will be the free parameters. Given the results of Paper II on the evolution 
of deep ISOCAM counts, it is expected that at z < 1.5 the starburst galaxies will have a stronger 
evolution rate than what has been assumed for normal late- type galaxies (assumption (1)) and for 
galaxy with AGNs (assumption (2)), consistent with the preliminary identifications of faint ISOCAM 
sources (Flores ct al. 1999). 

(4) Again for the sake of simplicity, we assume that all above mentioned evolution rates will have the same 
behavior after z — z\ 1 namely v\ = V2 = v% = q2- This is because beyond z~ 1.5 (the highest redshift 
detected for ISO galaxies), we know very little about the population of IR galaxies. 

4. Comparisons with Available Surveys and Constraints on Evolution Parameters 

4.1. Comparisons with ISOCAM 15/im Band Surveys: Constraints on Evolution of z < 1.5 

In what follows, we will compare our simulations with surveys in different bands found in the literature. 
The evolution models considered in this paper arc listed in Table 4. 

We start with the surveys in ISOCAM 15/um band (Elbaz et al. 1999; Serjeant et al. 2000), where the 
deepest and the most comprehensive ISO surveys have been conducted. Due to significant effects caused by 
the UIBs, certain features in the MIR counts can help to constrain the rate of luminosity evolution and that 
of density evolution separately (Paper I). If the sharp peak at /i5 Mm ~ OAmJy in the Euclidean normalized 
differential counts is indeed due to the UIBs in 6 — 8 /im (Paper II), which are redshifted into the 15/im 
band when z ~ 1, then a 15/im luminosity of ~ lO n L0 (yL nu at 15/im) could be inferred for a typical z=l 
ISOCAM galaxy. This imposes a strong constraint to the luminosity evolution rate of the major population 
of the ISOCAM sources, which under our assumptions (Section 3.3) are the starburst galaxies, at z ^, 1. 

In Fig. 7a and Fig. 7b, we compare the simulations of two evolution models, one has P2 = 2, u% = 4.2, 
Z\ = 1.5 and v\ = i>2 = i>3 = Q2 = —3 (Model 1, 'Peak Model') and the other p2 = 2, v,2 = 4.2, z\ = 1.5 and 
v\ = V2 = «3 = Q2 = (Model 2, 'Flat Model') with the number counts of ISOCAM 15/im data surveys 
(Elbaz et al. 1999; Serjeant et al. 2000). Both models fit the ISO data very well. Namely, although the 
'Flat Model' predicts about a factor of 2 more sources at /ls^m = O.OlmJy, there is little difference between 
the two models above the sensitivity limit of ISOCAM surveys ( 0.1m Jy). At 0.1m Jy the 15/zm counts are 
very insensitive to galaxy evolution at z > 1.5. 
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In Fig. 7c, 7d simulations of two other evolution models, Model 3 (p2 = 1, t*2 = 5, z\ = 1.5, 
«i = t>2 = U3 = q<i = —3) and Model 4 (p2 = 3, U2 = 3.5, £i = 1.5, Vi = i>2 = v% =92 = — 3) are also 
plotted. The former assumes more luminosity evolution (1L2 — 5) and less density evolution (j>2 = 1) for the 
starbursts at z < 1.5. It gives a slightly less good fit to the 15/um data (Fig. 7c), predicting a shallower peak 
at slightly brighter flux level than that of the data. The latter assumes less luminosity evolution (uz = 3.5) 
but more density evolution (p2 = 3) for the starbursts before z < 1.5. It fits the 15/zm data very well 
(Fig.7d). 

In Fig. 8, we compare the redshift distributions predicted by the 'Peak Model' model and by the 'Flat 
Model' with the data. Both the data of Aussel et al. (1999) and of Flores et al. (1999) are assumed to be 
complete at the 50% level (i.e. half of the sources are missing from the two plots due to lack of redshifts). 
By definition (Table 4), the 'Peak Model' and the 'Flat Model' differ only at z > 1.5. For the HDF-North 
survey, the model predicts a slightly higher median z (~ 1) compared to the median of the data (z ~ 0.7). 
However, since the redshift data are not complete and the high redshifts are more likely to be missing (more 
difficult to measure), and since the HDF-North is such a tiny field that any cluster at a given redshift (e.g. 
at z ~ 0.7) can affect the redshift distribution significantly, we feel this discrepancy is not inconsistent 
with our models. The CFRS survey is shallower (and wider) than the HDF-North survey, but its redshift 
distribution is more skewed toward the high z end, supporting our argument that the z distribution of the 
ISOCAM sources in the HDF-North field might not be representative for /ls^m > 0.1 sources. The 'Peak 
Model' predicts a median redshift of 0.75 for the CFRS survey, very close to that of the data. 

4.2. Comparisons with the CIB and the SCUBA 850/jto Band Surveys: Constraints on 

Evolution of z > 1.5 

Because of the negative K-correction in the submm bands, the best constraints on galaxy evolution 
beyond z ~ 2 come from the CIB and the submm counts. Our model predictions for the CIB are derived 
by summering up the flux densities of all sources in a very deep simulation (/24^ m > I0~ 9 mjy), for the 
band in question. In Fig. 9 we compare the CIB predicted by four different models to the observations. The 
solid curve is obtained by the 'Peak Model', which fits the CIB very well (except for the 60/j,m point of 
Finkbcincr et al. 2000). The dotted curve is the result of the 'Flat Model', which over predicts the submm 
CIB substantially (by a factor of ~ 2). The stronger density evolution model in Fig. 7c (Model 4, dot-dashed 
curve in Fig. 9) over-predicts the CIB by ~ 30 — 50%. In particular, it marginally violates the upper limits 
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set by the TeV gamma-ray observations (Stancv and Franceschini 1998) in the MIR. Another model (Model 
5) which is otherwise same as the 'Peak Model' except for a less steep drop (v\ = v 2 = v$ = q 2 = —1.5 
instead of v\ = v 2 = v 3 = q 2 = —3) after z=1.5, is also plotted in Fig. 9 (dashed curve). It slightly 
over-predicts the CIB around SOOfim. 

In Fig. 10 we compare results from simulations of the 'Peak Model' and of the 'Flat Model' to the 
850/j,m SCUBA counts. As in the case of the CIB comparisons, the 'Peak Model' fits the data very well, 
while the predictions of the 'Flat Model' are significantly higher than the data. 

4.3. Comparisons between Best-fit Model and the Surveys in Other Bands 

The agreements between the predictions of our best-fit model, the 'Peak Model', and the data from 
the ISO surveys at 90/im (Fig. 11) and at 175^m (Fig. 12) are very good. Note that near the faint ends of 
the 175/irn counts, incompleteness at level of £ 50% is expected (Dole et al. 2000). Corrections for this 
incompleteness in those data will make the agreement between our model predictions and the data even 
better. 

These good agreements may not be very surprising given that the model plotted here is mostly 
constrained by fitting the 15/im survey data (Fig. 7), and that the three FIR bands are linked to the 15/im 
band by a SED model that is very robust (Section 2). In Fig. 13 the predictions of our best-fit model are 
compared with the IRAS surveys at 60/im (Fig. 13). Here the data show a large spread, and our model is 
about 10% higher relative to the data at log fea < 0.3 (Jy) if the Grcgorich ct al. (1995) data are ignored 
(they may suffer from over-estimated; Bertin et al. (1997)). If real, the small over-prediction of these bright 
counts might be a consequence of a high local normalization (cf. Lonsdale et al. 1990). We will investigate 
this possibility further in our next paper. 

We compare the redshift distribution of IRAS 60/im sources observed by Oliver et al. (1996) and that 
predicted by our best-fit model in Fig. 14. The overall agreement looks quite good. On the other hand, 
there seems to be a slight excess of data points in the low z region (z < 0.014) and the opposite in high z 
(z > 0.033) region. It is not clear whether these are due to large scale structures (the amplitudes of the 
deviations are of the same order as the fluctuations due to the large scale structures) or they hint that the 
luminosity evolution rate of the best-fit model is too strong. It is also not clear whether the low fraction of 
high z galaxies could be explained, at least partially, by the incompleteness of the redshift survey (~ 10%, 
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Oliver et al. 1996). Future deeper surveys will certainly help to answer these questions. 

Finally we present the predicted contributions of IR-bright sources to the counts in the NIR K band, 
the optical B band, and the radio continuum 20cm band, and compare them with survey data from 
the literature. Again only the predictions of our best-fit model, the 'Peak Model', are considered here. 
Euclidean normalized differential counts in K band arc plotted in Fig. 15. Compared to the data points 
taken from various K band surveys (Soifer et al. 1994; Gardner et al. 1996; Bershady et al. 1998; Minezaki 
et al. 1998), the 'Peak Model' predicts that for K < 22, about 50 to 80% of the sources are IR bright, in 
reasonable agreement with the fraction of late-type galaxies among K band sources (Huang et al. 1998). It 
appears that the model prediction is significantly below the counts of K > 22, though there arc only two 
data points there. A tentative comparison with HST NICMOS H band deep counts (Storrie-Lombardi, 
private communication) shows that the 'Peak Model' predicts about 50% of the NIR counts down to H=26.5 
mag. 

In Fig. 16, we compare predicted contributions from IR sources to the B band counts. Since the E/SO 
galaxies constituent only ~ 20% of the optical galaxies (Glazebrook et al. 1995), which is a general result 
holding even for a very faint sample down to mi ~ 24.25 mag (Driver ct al. 1995), IR-bright sources should 
dominate the B band counts. The discrepancy near the bright end (B ~ 15) could be explained by the 
local density enhancement (local super-cluster). At the faint end (B > 25), the 'Peak Model' also slightly 
undcr-predicts the counts by ~ 50%. It should be noted that at these faint B band flux levels, where many 
sources have high redshifts and the B band flux is actually due to the rest frame UV emission, the model 
predictions suffer large uncertainties because the UV SEDs used in the calculation are not well constrained 
(Section 2). This will be the focus of our next paper. 

In the radio continuum 20 cm band, the bright sources (> lmJy) are mostly due to the early-type radio 
galaxies (Condon 1984) which are not IR emitters, while the faint, sub-mJy sources are mostly late-type 
galaxies (Condon 1984). This is indeed what the 'Peak Model' predicts (Fig. 17): at flux levels brighter than 
1 mJy, the IR sources contribute less than 10% of radio counts. At faint flux levels (~ O.lmJy), they can 
fully account for the radio counts. 

5. Predictions for Multi-Band SIRTF Surveys 
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5.1. Number Counts and Confusion limits 

In this section we will make predictions using our best-fit model ('Peak Model'). We will concentrate 
on future surveys with SIRTF, which will be launched in mid 2002. All three MIPS bands (24/im, 70/im, 
160/im) and all four IRAC bands (3.6/im, 4.5/im, 5.8/im and 8/im) bands are considered. For the sake of 
simplicity, we assume all SIRTF bandpasses are 10% (X/SX = 10)f]. We note that, since the E/SO galaxies 
are not included in the model, we may significantly underpredict the counts in the IRAC bands, particularly 
for the shorter wavelength bands (i.e. the 3.6/Ltm, 4.5/fm, and 5.8/xm bands). Therefore, our predictions for 
the counts and the confusion limits in these bands should be treated as lower limits. 

In Fig. 18 we plot the predicted integral counts for the three MIPS bands and three IRAC bands 
(3.6/im, 5.8/im and 8/im). In all MIPS bands, the counts are dominated by the starburst component 
(including any heavily obscured AGN). In the IRAC bands, normal spirals and galaxies with AGNs give 
significant contributions to counts brighter than 0.3mJy (in the 3.6/im and 5.8fim bands they out-number 
the starbursts). At fainter flux levels (< O.lmJy) the starburst component dominates the counts in the 
IRAC bands as well. 

Confusion limits shown in Fig. 18 have been computed from the predicted number counts, assuming 
parameters appropriate for SIRTF. The method used is the same as that used by Hacking and Soifer 
(1991)-namely, summing (in quadrature) the contributions from all sources fainter than the limit within 
the beam, via numerical integration of Eq. (19) in Hacking and Houck (1987). An Airy function computed 
from the nominal wavelength is used for the beam. For wavelengths greater than 6/xw, a telescope diameter 
of 85 cm is used in accordance with the required diffraction-limited performance of the SIRTF telescope. 
For wavelengths shorter than 6/xm, an effective telescope diameter is used to make the beamsize larger. 
Specifically, an effective diameter of 82 cm is used for 5.8/im, and 51 cm is used for 3.6/im. As noted above, 
the confusion limits for the short wavelength IRAC bands such as the 5.8/im and the 3.6/im bands should 
be treated as lower limits because of the omission of the contribution from counts due to E/SO galaxies. 
Also for the MIPS 160/im band, confusion caused by Galactic cirrus may be significant, especially in high 
cirrus regions (Gautier et al. 1992; Helou and Beichman 1990). 

In Fig. 19 redshift distributions predicted by the best-fit model ('Peak Model') and by the 'Flat 



3 See "Space Infrared Telescope Facility (SIRTF) — Observer's Manual (Version 1.0)" for details of SIRTF 
instruments. 
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Model' are plotted. For a 24/im survey limited at jiAjxm = 0.055mJy (5<7 con /), the best-fit model predicts 
a prominent peak at about z=l. The starburst component overwhelmingly dominates the counts at 
z > 0.5, while the normal galaxies dominate the small z end. The contribution from galaxies with AGNs 
is never important (< 10%), although it could be significantly higher for obscured AGNs classified here as 
"starbusrts" due to their colors. The shape of the redshift distribution predicted by the 'Flat Model' is very 
different, with a second peak at z ~ 2, caused mainly by the prominent UIB features between 6 — 8/im 
(Paper I). This peak is cut off in the 'Peak Model' because of the steep decrease of the star- formation rate 
at z>1.5. However, the turn-over redshift z\, which has been assumed to be 1.5 in all our models presented 
here, is poorly constrained because ISO surveys are too shallow whereas the SCUBA and CIB are more 
sensitive to sources at z > 3. If in reality z\ > 2, then it will be revealed by the second peak at z ~ 2 in the 
redshift distribution of SIRTF 24/im sources. 

5.2. SIRTF Color-Color Diagrams 

In Fig. 20a and Fig. 20b we plot the ji^ml fi\xm vs. /70/jm//24^m color-color diagrams predicted by 
the 'Peak Model' and the 'Flat Model', respectively. Simulations of both models are confined to 1 deg 2 
with the following flux limits: ,/*24/im > 0.055'mJr/ (5<T c „ n /) and fro^m > 2.6mJy (2er con /). Some general 
trends are visible in these plots. Galaxies with AGNs are mostly in the relatively low /70/xm //24/im ratio 
region (log(/ 70AJm / / 2 4 M m) ~ 1)), while starburst galaxies dominate the high / 7 o-im//24/im ratio region 
(log(/7o Mm //24^m) ~ !))• Normal galaxies, with relatively low luminosities and therefore seen only with 
z < 1, are concentrated in the low /24/im/ fa^m ratio region. It should be pointed out that these trends 
are closely related to our definitions of the three populations (Section 2), which are separated according to 
the IRAS color ratios of feoum/ j '25 nm an d /25/jm//i2/im- While the criteria adopted are based on empirical 
correlations (Fang et al. 1998; Paper I), the clear boundaries in the definitions of the different populations 
will have artificially enhanced the trend shown here. 

In Fig. 21a and Fig. 21b we plot the predictions by the same two models for the MIPS color-color 
diagram, namely /i60/im/ /70/im vs. /70/im/ /24/im- The symbols are the same as in Fig. 20. In these 
plots, in addition to the flux cut-offs of /24/tm > 0.055mJy and fro^m > 2.6mJy, it is also required that 
/i60/oti > 38m Jy (2a con f). The simulations cover 10 deg 2 . Most high z galaxies (large crosses in the plots) 
are concentrated in the high fwo^m/f7ap, m end (log(/i 6 o^m//70^m) ~ 0.7), due to the K-correction effect 
associated with the curvature in the 20 — 160/xto wavelength range in most of the SEDs in our library. 
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The apparent 'tracks' in these figures are a result of the manner in which an SED is selected for each 
simulated galaxy from the SED library. As luminosity increases with increasing z (for luminosity evolution 
models) a model galaxy is assigned a random SED from the relevant luminosity bin. Some bins contain 
small numbers of SEDs so the 'tracks' reflect single SEDs randomly assigned to model galaxies at higher 
and higher z (and therefore L). The 'tracks' can discontinue when model populations 'jump' to higher SED 
library luminosity bins. Although small random scatter could be incorporated into these model colors, we 
have chosen not to do that so that the figures can be more easily analyzed. 

The IRAC color-color diagrams, such as /ks/W/Wm vs. fs fJ , m /h.s f j,m diagram, are affected by the 
exclusion of E/SO galaxies in our current models. Also the NIR spectral features considered by Simpson and 
Eisenhardt (1999), e.g. the 1.6/ito H opacity minimum and the 2.3/ito CO bandhead, are not considered 
in our SED model (Section 2), therefore the changes in the IRAC colors due to the K-correction caused by 
these features are absent in our model predictions. Improvements in these aspects will be addressed in our 
next paper. 

6. Discussion 
6.1. Comparisons with Previous Models in Literature 

6.1.1. Comparisons with Models in Paper I and Paper II 

Three new features separate the models presented here from the previous models in Paper I and Paper 
II: 

(1) An approach even more empirical than that used in Paper I is adopted in this work. Realistic SEDs 

of local galaxies are attached to sources of all redshifts. This not only enabled reliable K-corrections, 
but also linked the surveys in different bands coherently. Because of this feature, our model is the 
first in the literature that can predict the correlations and the dispersions in color-color diagrams for 
multi-band surveys. 

(2) The MIR-FIR SEDs in the SED sample are extended to a much wider range covering from UV (1000A) 

to radio (20cm). This links the model predictions for IR surveys to the optical and radio continuum 
surveys. Because of the wide wavelength coverage, our model can calculate the contributions from the 
IR sources to the cosmic background from optical to radio bands. 
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(3) Instead of evolving all IR sources as a single population, they are divided into three populations 
(normal late- types, starbursts, and galaxies with AGN) which are assumed to evolve differently. Given 
the observational evidence (Aussel et al. 1999; Elbaz ct al. 1999; Flores et al. 1999), most of the 
evolution is attributed to the starburst (interacting) galaxies. 

In Paper II, it was found that a pure luminosity evolution of form L oc (1 + z) 4 ' 5 (z < 1.5) can fit 
the ISOCAM 15pm surveys and IRAS 60pm surveys best (both surveys are insensitive to the evolution 
at z > 1.5). This evolution rate is significantly stronger than those derived in previous IRAS studies 
(L oc (1 + z) 3 ), while at same time it confirms that luminosity evolution models fit the data better than 
density evolution models (e.g. Pearson and Rowan- Robinson 1996). 

In this paper, while assuming that normal late-types and galaxies with AGNs undergo pure luminosity 
evolution (L oc (1 + z) and L oc (1 + z) 3 ' 5 , respectively), we found that a model assuming that both the 
luminosity and the density of starburst galaxies evolve significantly [L oc (1 + z) 4 ' 2 and p oc (1 + z) 2 ) before 
the turn-over redshift (z=1.5), fits the data best. A model with more luminosity evolution and less density 
evolution (L oc (1 + z) 5 and p oc (1 + z)) predicted a peak in the Euclidean normalized differential counts 
at 15/zm too shallow and too bright compared to the data (fig. 7c). Another model with less luminosity 
evolution and more density evolution (L oc (1 + z) 3,5 and p oc (1 + z) 3 ), while fitting well the 15pm counts, 
predicted too much IR background (fig.7d and Fig. 8a). 

The prediction of our best-fit model of a density evolution of starburst galaxies on the order of (1 + z) 2 
is consistent with the theoretical prediction (Carlbcrg ct al. 1994) and with previous observations (Infante 
et al. 1996; Roche and Eales 1999; Le Fevre et al. 2000) of the evolution of merger rate. It is also very 
consistent with the close relation between starbursts and galaxy interactions/mergers (see Sanders and 
Mirabel 1996 for a review). 

The constraints on the evolution of IR sources at z > 1.5 are set by the CIB and the 850pm SCUBA 
counts. Because of the negative K-correction, high redshift galaxies contribute significantly to the source 
counts and the cosmic background in the submm bands. It is found that in order to avoid over-predicting 
the CIB and the SCUBA counts, a steep decrease after z=1.5, in the form of (1 + z)~ 3 , in the evolution of 
IR sources is required. This is different from the result in Paper II, where the 'Flat Model' (luminosity and 
density remain constant after z=1.5) can fit the CIB satisfactorily (Fig. 7 of Paper II). This is due to the 
difference in the submm SEDs adopted in these two papers. In Paper II a single template for the submm 
SEDs is adopted for all sources, specified by T=40K and (3 — 1.5 (Blain et al. 1999). In this paper, the 
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submm SEDs are determined using empirical correlations between IRAS fluxes and the submm fluxes. This 
results in a higher submm luminosity for a given MIR luminosity (Li^m)- Since little is known about the 
IR-submm SEDs of high z sources, and we cannot test our assumption that the correlations between the 
IR and submm fluxes of local galaxies can be applied to high z galaxies, the rapid decrease of the SFR 
predicted by our best-fit model is highly uncertain. 

Because of the steep decrease after z=1.5, our best-fit model predicts a redshift distribution with a 
single peak at z^ 1 for a survey with f2ifj.n1 = O.ObbmJy (5<7 C(m /). This is very different from the double 
peaked distribution predicted by the 'Flat Model' (Fig. 19, see also Fig. 14c of Paper I). This provides a 
straightforward test to distinguish these models once the redshifts of a sample of SIRTF 24/im sources are 
obtained. 

6.1.2. Comparisons with Roche and Eales (1999) 

In the multi-IR-band model of Roche and Eales (1999), a density evolution of p oc (1 + z) 2 is also 
assumed. However, the luminosity evolutionary rate of starburst galaxies found by Roche and Eales (1999), 
L oc (1 + z) 2 (z < 1), is significantly weaker than predicted by our best-fit model: L oc (1 + z) 4 - 2 {z < 1.5). 

6.1.3. Comparisons with Dole et al. (2000) 

Compared to the other two populations (normal late- types and galaxies with AGNs), our best-fit model 
predicts much stronger evolution for z<1.5 in the starburst galaxies population. This is similar to the 
model of Dole et al. (2000). However, they confine their starburst component, called 'ULIRGs', to galaxies 
with > 2 x 10 11 Lq. In order to fit the ISO counts they have to evolve the 'ULIRG' population significantly, 
resulting in a 'break' in the evolved luminosity function at z=2.5 (Fig. 3a of Dole et al. 2000). In Fig. 22, the 
evolved luminosity functions predicted by our best-fit model are plotted. At z>l, the starburst component 
dominates the LF in the whole luminosity range (10 7 — 10 12 ). Our model does not predict a break such as 
presented in the model of Dole et al. (2000) . 
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6.1. 4- Comparisons with Blain et al. (1999) 

Blain ct al. (1999) used IRAS 60/ito counts, early results of ISO 175/ito counts, SCUBA 850/im 
counts, and the CIB to constrain several families of models for the evolution of IR sources. Since different 
cosmologies are used (fl m — 0.3, fl\ = 0.7 in this work, Cl m — 1, J7a = in Blain ct al. 1999), it is difficult 
to make a quantitative comparison. In the three panels in Fig. 23, we plot the cosmic luminosity density 
evolution from our best-fit model. In the upper panel, the MIR luminosity density, calculated using the 
parameters of the best-fit model, is plotted versus the time since the Big-Bang (in fraction of the age of the 
Universe to). In the middle panel, the density of bolometric luminosity (0.1 — 1000 fim) obtained from a 
simulation (,/24/jm > 10~ 9 mJy, 0.01 deg 2 ) based on the best-fit model is plotted. In the lower panel, results 
for the density of IR luminosity (3 — 1000 fim) from the same simulation is plotted. For IR sources, the IR 
luminosity density is proportional to the star formation rate per unit cosmic comoving volume. 

Comparing the lower panel of Fig. 23 with Fig. 9 of Blain et al., our best-fit model is closer to their 
'Peak-Models' (rise — peak — drop) than the 'Anvil-Models' (rise — peak — flat). As pointed out by Blain et 
al. (1999), if indeed the contributions from AGNs to the CIB and to the submm counts are negligible, the 
'Peak-Models' are favored because the 'Anvil-Models' may over predict the metal contents of the Universe. 

6.1.5. Comparisons with Rowan-Robinson (2001) 

Recently Rowan- Robinson (2001, hereafter RR01) developed a model to study the evolution of galaxies 
using multi-band IR surveys. This goal is similar to that which motivated this paper. Consequently, there 
are many similarities between RR01 and the work presented here. 

(1) Strategy: Both RR01 and we use all the available surveys to constrain the evolution, then use the 

best-fit model to make predictions for future surveys. 

(2) Methodology: Both RR01 and we use local SEDs to link all surveys together, though RR01 started 

from the IRAS 60/ito band and we started from IRAS 25/ito band. 

(3) Results: Both the models of RR01 and of this paper can fit all ISO surveys, SCUBA 850/ito surveys, 

and the CIB. 

(4) These agreements are the natural consequence of the fact that in both papers the model parameters 

are tuned to fit the counts in these bands. In Fig. 24a, a comparison between the model predictions on 
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the cosmic luminosity density evolution by RR01 (the model for A = 0.7, dot-dot-dot-dashed curve) 
and by our best-fit model ('Peak Model', solid curve) are plotted. Though the function forms are 
different, the overall trends arc quite similar. 

On the other hand, there are some major differences between RR01 and this paper: 

1) Assumptions: RR01 assumed that all the IR sources in his model, which he decomposed into four 

spectral components (cirrus, M82-like starburst, Arp220-likc starburst, and AGN dust torus), evolve 
as a single population. We have assumed that the three populations (normal late-types, starbursts, 
and galaxies with AGNs) in our model evolve differently. RR01 also adopted a different functional 
form (exponential) to describe the evolution, while the evolution functions in this paper are simple 
power-laws. 

2) SEDs: We have used 800+ empirically determined SEDs, which arc binned according to population and 

luminosity, to constrain the K-corrections and the links between counts in different bands. This is 
more sophisticated than the approach of RR01, where four theoretical SEDs taken from the literature 
(Efstathiou et al. 2000b; Yoshii and Takahara 1988) are adopted for the four populations, respectively. 

3) Evolution of SEDs: In RR01 models, the SEDs of IR sources undergoing luminosity evolution do not 

evolve. This is different from our assumption for the luminosity dependence of SEDs, namely when 
the luminosity of the IR sources evolves, the SED changes with the luminosity (see Section 6.3). 

4) Predictions for redshift distribution of SIRTF sources: In Fig. 24 we compare the redshift distribution 

of sources in the SIRTF 70/im band (/70/jm > 5 mJy) predicted by RR01 and by our best-fit model. 
While RR01 predict that most of these sources ('cirrus' and starbursts) have z < 0.5 (median 
~ 0.4), our best-fit model predicts the opposite, namely most of sources (starbursts) have z > 0.5 
(median <~ 1). This difference is mainly due to the different evolution rates of normal-spirals ('cirrus 
galaxies') in the two models: In our best-fit model, the evolution rate of the normals is rather low 
(evolution index = 1.5). Since the normals dominate sources in the local Universe, the evolution 
rate of overall IR sources at low z is low in our best-fit model. Only when the starbursts overtake 
the normals as the dominant population, which occurs at z ~ 0.5, does the high evolution rate 
of the starbursts start to significantly affect the IR sources on the whole. On the other hand, in 
RR01, all populations evolve with the same rate. So the evolution index is ~ 3 at small redshifts for 
normal spirals. Consequently, as shown in the upper panel of Fig. 23, the star formation rate (roughly 
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proportional to the luminosity density) of RR01 is about a factor of 2 higher than our best-fit model 
at z < 0.5. The SFR predicted by our best-fit model starts to catch up with RR01 at larger z and 
then overshoots until z~ 1.5. It follows that a redshift survey of SIRTF sources in the 70/im band will 
serve as a good test to distinguish the models. It should be noted that for the sample of faint IRAS 
60/iTO sources (feoum > 200 mJy, Oliver et al. 1996), our best- fit model predicts too many (about 
a factor of 10) sources with z> 0.4 compared to the observations (6 out of 1400: Rowan-Robinson, 
private communication). However, it is not very clear to what extent the discrepancy is due to biases 
introduced by large scale structures and by the incompleteness of the redshift survey (see the end of 
Section 4). RR01, which predicts about a factor of 2 less such sources than our best-fit model, gives a 
better fit to the redshift data. 

5) CIB: The difference in the evolution rates in different populations between RR01 and our best-fit 
model results in a significant discrepancy on the relative contributions from different populations 
to the counts and to the CIB. As shown in Fig. 25, our model predicts that the IR-submm CIB is 
predominantly due to the starburst galaxies. The best-fit model of RR01 predicts that most of the 
CIB is due to 'cirrus galaxies' (close to 'normals' in our model) which, because they are assumed to 
evolve as strongly as starburst galaxies, maintain their dominance to the IR emission from the local 
Universe throughout the high-z Universe. Both of our models predict negligible contributions from 
AGNs to the CIB (< 10% at any given wavelength), although the 'starburst' model populations could 
include a significant heavily obscured segment of the AGN population that are thought to contribute 
significantly to the Cosmic Xray Background (XRB; eg. Gilli, Salvati and Hasinger 2001). 

6.2. Star Formation History 

Much attention has been attracted since Madau et al. (1996) related source counts and reshift 
distributions obtained from deep UV/optical surveys to the star-formation/metal-production history of 
the Universe. Since then the so-called "Madau-Diagram" has been revised many times through various 
improvements, including (1) the corrections for the effect of dust extinction on UV/optical luminosities 
(Madau et al. 1998; Pozzetti et al. 1998, Steidel et al. 1999), (2) results from less extinction sensitive 
Balmer line surveys (Gallego et al. 1995; Tresse and Maddox 1998; Yan et al. 1999; Glazebrook et al. 
1999), (3) results from submm SCUBA surveys (Barger et al. 1999). 

Adopting the conversion factor of Kennicutt (1999), SFR (M Q yr _1 )=L/# x 4.5 10~ 44 (erg s _1 ), we 
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convert the IR (3 — lOOO^tm) luminosity density curve plotted in Fig. 23 (lower panel) to a star formation 
rate curve. In Fig. 26 the result is compared with the survey data found in the literature. Not surprisingly, 
the model prediction is in very good agreement with the results from ISOCAM surveys (Flores et al. 1999) 
because at z < 1.5 the model is constrained by the ISOCAM data. On the other hand, the model prediction 
is slightly higher than the results of the UV/optical and Ha surveys of z < 1.5. After the turn-over 
redshift, our best-fit model predicts a quick decrease, similar to what have been suggested by the UV data 
of Lyman-break galaxies. 

The 'Flat Model', represented by the dashed line, predicts too much CIB (Fig. 9) and too many counts 
in the SCUBA 850^m band (Fig. 10). The SFR of z > 3 galaxies predicted by the 'Flat Model' is about a 
factor of 7 higher than those derived from the UV data of Lyman-break galaxies, even after the UV data 
are corrected for the extinction (about a factor of 3, Stcidel et al. 1999). Although we can not rule out 
the 'Flat Model' because of the uncertainties associated with the evolution of the submm SEDs, and those 
associated with the extinction corrections of the UV data, it is certainly disfavored by our results. On the 
other hand models between the 'Peak Model' and the 'Flat Model', such as the Model 5 plotted by the 
dotted line (see also Fig. 9), are certainly allowed by our results. 

6.3. Uncertainties Introduced by Model Assumptions 

The most important assumption in this work concerns the applicability of the SED vs. IR luminosity 
relation, found for local IR galaxies (Soifer and Neugebauer 1991; Fang et al. 1998), for high z galaxies. 
Namely we assume that the SEDs do not evolve with time for a given luminosity. Note that this is different 
from assuming that the SEDs of galaxies do not evolve at all, because our models indeed predict strong 
luminosity evolution for IR galaxies and therefore, under our assumption, more galaxies in the early epochs 
have SEDs similar to local luminous IR galaxies which are significantly different from the SEDs of local L* 
galaxies. 

This assumption is not very well constrained. Not much is known about the IR SEDs of high z galaxies, 
in particular no current FIR (A ^ 30/xto) instrument is sensitive enough to detect them in the FIR bands. A 
few z ^ 1 galaxies have been identified in the ISOCAM surveys (predominantly in ISOCAM 15/im surveys, 
Aussel et al. et al. 1999; Flores et al. 1999), and the optical SEDs of some of them have been obtained 
using HST images (Flores et al. 1999). The submm SCUBA surveys (Hughes et al. 1998; Barger et al. 
1998; Blain et al. 1999) may have detected quite a few high z galaxies. But the positional uncertainties of 
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the SCUBA sources are so large that it is very difficult to make cross identifications in other bands. In fact 
so far only a handful (3 - 5, Frayer: private communication) high z SCUBA sources have published high 
confidence optical identifications and redshifts (Ivison et al. 1998; Barger et al. 1999; Frayer et al. 1999a). 
The most recent results of Ivison et al. (2000) show that indeed the SEDs of these galaxies have similar 
shapes to local ULIRGs such as Arp 220 or Mrk 231. 

High-z galaxies seen in deep surveys are high luminosity galaxies. Thus, to the extent that they have 
similar SEDs to their local counterparts, our results will be valid. Namely, intrinsically faint galaxies at 
high z may have different SEDs than their local counterparts, but this won't have any significant effect on 
the predictions for number counts and the CIB. It is likely that in high z ULIRGs the luminosities in all 
bands are predominantly radiated in localized starburst regions and/or AGNs, where the physical processes 
determining the luminosity and the SED (e.g. AGN related processes, star formation, radiative transfer, 
etc.) are similar to those in local ULIRGs. In particular, both the theoretical arguments and observational 
evidence (Soifer et al. 1998; Armus et al. 1998) show that dust is produced rapidly after the first star 
formation episode in galaxies. Therefore those high luminosity galactic nuclei in the early Universe are 
likely to be optically thick (Soifer et al. 1998; Armus et al. 1998), similar to their local counterparts. 

It is expected that for large samples of galaxies SEDs covering 4000A — 160/zm will be available when 
SIRTF deep surveys and follow-up optical/NIR surveys are carried out. Then the assumption that the SED 
vs. IR luminosity relation of local IR galaxies also holds for high z galaxies will be fully tested. 

Throughout this paper, we have assumed the A cosmology (£1a = 0.7, Q m = 0.3), which is favored 
by recent observations of type I supernovae in distant galaxies (Perlmutter et al. 1997; Garnavich et al. 
1998). Compared to the standard Q\ = 0, Q m = 1 Einstein-de Sitter cosmology, the co-moving volume 
corresponding to a given z is significantly larger in the A cosmology (Carroll et al. 1992). Althogh this 
volume factor is partially balanced by relatively fainter flux for a given luminosity and given z in the A 
cosmology, it still affects significantly the predictions of the SCUBA counts and the submm CIB, where the 
contribution from high z sources is large. This is one of the reasons why a steep decrease after z=1.5 is 
favored by our best-fit model (Fig. 26). If the Einstein-de Sitter cosmology were adopted, a flatter SFR at 
early epochs would have been obtained. 

The assumptions for the evolution rates of normal late-type galaxies and of galaxies with AGNs 
(Section 3.3) are based on the results of optical and NIR surveys of these sources. Whether they also apply 
to the IR bands will have to be tested with future IR surveys (e.g. SIRTF surveys). Nevertheless, if indeed 
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most of the evolution in IR bright galaxies is due to starburst galaxies (and possibly also highly obscured 
AGN), as suggested by ISOCAM surveys (Elbaz et al. 1999; Flores et al. 1999), these uncertainties will 
have minimal effects on our results for the overall IR evolution. 

Finally, the prominent peak at z = 1.5 in the SFR v.s. redshift plot (Fig. 25) is an artifact due to the 
simple two-step power-law function form adopted for the evolution functions (Eq(7) and Eq(8)). However 
the SFR at z<1.5 is mostly constrained by the ISOCAM 15/ito surveys, which indeed reach as deep as 
z=1.5. As was argued in Paper II, a very strong evolution all the way back to z~1.5 is needed to explain 
the sharp peak at /i5 Mm ~ 0.4 mJy in the Euclidean normalized differential 15/im counts (Elbaz 1999). 

7. Conclusions 

We have developed empirical 'backward-evolution' models for multiband IR surveys. A new Monte- 
Carlo algorithm is developed for this task. It exploits a large library consisting of realistic Spectral Energy 
Distributions (SEDs) of 837 local IR galaxies (IRAS 25/zm selected) from the UV (1000A) to the radio 
(20cm), including ISO-measured 3-13/im unidentified broad features (UIBs). The basic assumption is that 
the local correlation between SEDs and Mid-Infrared (MIR) luminosities can be applied to earlier epochs 
of the Universe. A A-Cosmology (fl\ = 0.7, Cl m = 0.3) has been assumed through out the paper. By 
attaching an SED appropriately drawn from the SED library to every source predicted by a given model, 
the algorithm enables simultaneous comparisons with multiple surveys in a wide range of wavebands. 

IR galaxies are divided into three populations: (1) normal late type galaxies ('normals'), (2) 
starburst/interacting galaxies ('starbursts'), and galaxies with AGNs ('AGNs'). Different cosmic evolution 
is assumed for these different populations. Parameterized (power-law) luminosity evolution functions 
(Fi(z)) and density evolution functions (Gi{z)) of the form 

Fi(z) = (l + z) Ui (z< Zl ) 

= (1 + z) Vi {z 1<Z < z ) 

Gi(z) = (l + z) Pi (z< Zl ) 

= (1 + z) qi ( Zl <z< z ) 

are adopted, with z\ — 1.5 and zq — 7. At z < 1.5, for 'normals' (i=l) and 'AGNs' (i=3) it is assumed 
Mi = 1.5 pi = (passive luminosity evolution) and w.3 = 3.5 yz = (evolution rate of optical QSOs), 



28 



respectively. The evolution rate of 'starbursts' (i=2) at z < 1.5 is determined by fitting the ISOCAM 15/um 
surveys. The best fit results are U2 = 4.2 p2 = 2. At 1.5 < z < 7, the evolution of 1R galaxies is mostly 
constrained by the submm counts and the CfB. ft is found that a 'Peak Model', with the pi and Ui values 
described above and v\ = v^ = V3 = qi = —3, gives the best fit. The 'Flat Model', which is the same as the 
'Peak Model' at z < 1.5 but v\ — v<i = V3 = qi = at 1.5 < z < 7, over-predicts significantly the SCUBA 
counts and the CIB. 

Remarkably, the best-fit model ('Peak Model') gives good fits simultaneously to all data (both number 
counts and redshift distributions) obtained from IR surveys, including ISOCAM 15/um, fSOPHOT 90/im, 
175/irn, IRAS 60iim, and SCUBA 850/ito. Predictions for contributions of IR-bright sources to counts in 
other wavebands, such as the optical B band, the NIR K band, and the radio continuum 20cm band, are 
also in agreement with the literature. This suggests that the model is robust. 

Predictions for number counts, confusion limits, redshift distributions, and color-color diagrams are 
made for multiband surveys using the upcoming S1RTF satellite. It is found that several S1RTF colors can 
be useful indicators of galaxy populations and redshifts. 



Appendix 



A. Correlations between Flux Densities in the FIR and Submm Bands 

In this appendix, empirical correlations between flux densities in the IRAS 60/xrn and 100 [im bands 
and in the submm bands at 170/im, 2A0/im, 450/im, 850/im, and 1200/zm are studied. The results have 
been applied to the SED model in order to extend the SEDs to the submm (up to 1200 /im) wavebands 
(Section 2). 

Flux densities in the submm wavebands were collected from the literature. In the 170/im band 
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(including all data between the 160^m and 180/^m bands), data for 29 galaxies were found in Devereux 
and Young (1992), Engargiola (1991), Klaas et al. (1997), Hippelein ct al. (1996), and Stickel et al. (1998). 
In the 240/zm band, flux densities were found for 20 galaxies (including 14 upper limits) in Odcnwald et 
al. (1998). In the 450fj,m band, 16 detections and 7 upper limits were found in Rigopoulou ct al. (1996), 
Frayer et al. (1999b), Eales et al. (1989), Dunne et al. (2000), and Alton et al. (1998). In the 850/^m 
band (including also 800 /im band), 121 flux densities (including 1 upper limit) were found in Lisenfeld et 
al. (2000), Dunne ct al. (2000), Alton ct al. (1998), Rigopoulou ct al. (1996), Hughes et al. (1990), Frayer 
et al., (1999b), Eales et al. (1989). In the 1200/im band (including all observations between the 1000/im 
and 1300^m bands), 62 detections and 4 upper limits were found in Andreani and Franccschini (1996), 
Rigopoulou et al. (1996), Hughes et al. (1990), Eales et al. (1989), Chini et al. (1986), Fich and Hodge 
(1993). Color-color diagrams of log(/ s „6 mm ,j//ioo^m) versus log(/ioo^m/ feo^m), where f subm m,i is one of 
the submm flux densities listed above, are plotted in Fig. 27 — Fig. 32. In Fig. 31 and Fig. 32, color-color 
diagrams of log(/i 2 oo/wn//ioo^m) versus log(/ioo /im //60/ i m) are plotted for the same sources, with the 
sources in the list of Andreani and Franccschini (1996) (19 sources) being aperture corrected in two different 
ways as given in Andreani and Franccschini (1996), respectively. Linear relations in the form of 

l°g(/sti6mm,i/ ' f 100 fim) = Ai + Bi X log(/l00/^m//60jUm) (Al) 

were derived (eyeball) from these color-color diagrams. In Table Al, the resulting Ai and Bi are listed. For 
a given source (with detected IRAS flux densities f^o^m an d Aoo^m) hi the SED sample, the submm flux 
density in any of the given bands is then estimated using the relation: 

fog(/su6mm,i) = fog(/lOOMm) + A t + Bi X log(/ioo M m//60^m). (A2) 

It should be noted that the correlations plotted in the FIR/submm color-color diagrams are generally 
rather poor. There are often only a small number of sources in a plot, and many of them are uppcrlimits. 
The data are also very heterogeneous, obtained from observations with widely different apertures. 
Aperture corrections were included only when available in the literature. When there are more than one 
observations for a given source in a given band, the data from the observation with the largest aperture is 
taken. Consequently, the lines in the color-color plots which represent our best (eyeball) estimates of the 
fog(/stt&mm,i//ioo^m) versus log(/ioo^m/JWm) relations are very uncertain. This reflects the true situation 
in the current literature of extragalactic submm sources. Nevertheless, it appears that the submm SEDs 
derived using flux densities predicted by these relations agree well with the observations for a wide variety 
of extragalactic sources (Fig. 2). And the good agreement between the simulated and measured 850/xrn and 
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1.2mm luminosity functions (Fig. 5 and Fig. 6) further support the validity of these empirical relations. We 
have deliberately avoided any modelling in deriving the relations, keeping them purely empirical. They will 
set constraints to theoretical models of dust heating in galaxies (e.g. Silva et al. 1998; Dale et al. 2001; 
Popescuet al. 2000). 
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Table 1. Bins of SED library 



population 


t i25^m bin 


number of 




log(Li) 


log(L 2 ) 


sources 


Normals 


6 


8 


3 




8 


9 


81 




9 


9.4 


65 




9.4 


9.8 


63 




9.8 


10.2 


36 




10.2 


11 


18 


Starbursts 


6 


8 


2 




8 


9 


31 




9 


9.4 


46 




9.4 


9.8 


85 




9.8 


10.2 


80 




10.2 


10.6 


78 




10.6 


11 


41 




11 


12 


16 


AGNs 


6 


10 


63 




10 


12 


129 



^ L25/jm = z/L„(25/xm), in units of L Q . Luminosity bins are defined by log(Li) < log(L25 Mm ) < log(L2). 
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Table 3. Parameters of LLFs of Different Populations 



Population 


number 


a 


a 


L*/Lq 


C (normalization) 


Normals 


456 


0.482 


3.876 


5.7 x 10 9 


0.00035 


Starbursts 


643 


0.268 


2.230 


7.9 x 10 9 


0.00066 


AGNs 


356 


0.336 


1.691 


6.9 x 10 9 


0.000090 
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Table 4. Evolution models. Parameters zq and z\ specify the galaxy formation time and the 'turn-over' 
redshift, respectively; Ui and Vi (i — 1,2,3) specify the luminosity evolution functions (Eq(7)), and pi and 
qi (i = 1,2,3) specify the density evolution functions (Eq(8)). 

Normals Starbursts AGNs 



model zi z U\ vi p\ q\ u 2 v 2 p 2 q_i u 3 v 3 p 3 q 3 

Model 1+ 1.5 7 1.5 -3 4.2 -3 2 -3 3.5 -3 

Model 2* 1.5 7 1.5 4.2 2 3.5 

Model 3 1.5 7 1.5 -3 5-3 1 -3 3.5 -3 

Model 4 1.5 7 1.5 -3 3.5 -3 3 -3 3.5 -3 

Model 5 1.5 7 1.5 -1.5 4.2 -1.5 2 -1.5 3.5 -1.5 

t 'Peak Model'. 
* 'Flat Model'. 
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Table Al. Parameters of log{f S ubmm,i/fioo^ m ) vs- log(/ioo M m/ TWm) relation 



waveband 


detections 


upper limits 


A % 


B t 


Figure 


170/ito 


29 




-0.3 


1.0 


Fig.27 


240^m 


6 


14 


-0.7 


1.36 


Fig.28 


450/ito 


16 


7 


-1.58 


1.44 


Fig.29 


850^to 


121 


1 


-2.40 


1.45 


Fig.30 


1200/xto 


62 


4 


-2.90 


1.45 


Fig.31, Fig.32 
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Table 2. Average SEDs (I og ( f u / f 25^)) of IR sources in log(_L25 /U m) bins. One-cr dispersions are also given. 









Normal Late-types 






AGNs 


log(A) 






log(i 2 5 M m 


/ ' Lq) bins 






log(i 2 5 M m 


1 Lq) bins 


(jim) 


7±1 


8.5 ±0.5 


9.2 ±0.2 


9.6 ±0.2 


10.0 ±0.2 


10.6 ±0.4 


8±2 


11±1 


-1.00 


-1.31±0.35 


-1.78±0.57 


-2.20±0.48 


-2.31±0.52 


-2.65±0.42 


-2.80±0.45 


-2.68±0.67 


-3.30±0.57 


-0.76 


-1.05±0.44 


-1.57±0.44 


-1.94±0.44 


-2.05±0.45 


-2.39±0.41 


-2.49±0.44 


-2.42±0.69 


-3.02±0.57 


-0.51 


-0.82±0.41 


-1.21±0.37 


-1.56±0.36 


-1.68±0.34 


-1.99±0.33 


-2.11±0.35 


-2.07±0.55 


-2.63±0.48 


-0.27 


-0.58±0.33 


-0.80±0.35 


-1.12±0.29 


-1.27±0.24 


-1.54±0.25 


-1.69±0.25 


-1.68±0.40 


-2.21±0.40 


-0.03 


-0.35±0.28 


-0.42±0.39 


-0.72±0.30 


-0.88±0.22 


-1.11±0.22 


-1.29±0.20 


-1.31±0.38 


-1.81±0.39 


0.21 


-0.15±0.31 


-0.18±0.39 


-0.45±0.30 


-0.62±0.22 


-0.83±0.22 


-1.00±0.18 


-1.06±0.41 


-1.52±0.40 


0.46 


-0.42±0.29 


-0.47±0.31 


-0.67±0.23 


-0.79±0.18 


-0.96±0.18 


-1.08±0.14 


-1.19±0.31 


-1.42±0.34 


0.64 


-0.64±0.29 


-0.79±0.18 


-0.85±0.17 


-0.91±0.12 


-1.00±0.15 


-1.06±0.14 


-1.18±0.19 


-1.19±0.30 


0.83 


-0.27±0.29 


-0.43±0.18 


-0.49±0.17 


-0.57±0.13 


-0.69±0.15 


-0.75±0.14 


-0.89±0.19 


-0.93±0.27 


0.87 


0.16±0.28 


-0.02±0.18 


-0.10±0.17 


-0.19±0.13 


-0.37±0.18 


-0.54±0.15 


-0.69±0.25 


-0.83±0.25 


0.91 


0.05±0.29 


-0.10±0.18 


-0.16±0.16 


-0.22±0.12 


-0.28±0.14 


-0.25±0.14 


-0.64±0.21 


-0.70±0.26 


0.95 


-0.35±0.30 


-0.46±0.17 


-0.49±0.16 


-0.51±0.11 


-0.54±0.14 


-0.49±0.14 


-0.77±0.17 


-0.73±0.26 


0.99 


-0.43±0.29 


-0.59±0.17 


-0.65±0.16 


-0.69±0.11 


-0.74±0.13 


-0.74±0.12 


-0.80±0.18 


-0.78±0.26 


1.04 


-0.21±0.28 


-0.39±0.18 


-0.47±0.16 


-0.54±0.12 


-0.62±0.13 


-0.63±0.12 


-0.65±0.17 


-0.66±0.25 


1.08 


0.15±0.29 


0.00±0.17 


-0.06±0.16 


-0.12±0.11 


-0.19±0.14 


-0.17±0.13 


-0.43±0.17 


-0.48±0.24 


1.12 


0.07±0.29 


-0.07±0.17 


-0.11±0.15 


-0.15±0.11 


-0.19±0.13 


-0.18±0.11 


-0.39±0.15 


-0.40±0.22 


1.16 


-0.07±0.29 


-0.21±0.16 


-0.26±0.14 


-0.32±0.10 


-0.39±0.12 


-0.42±0.08 


-0.39±0.14 


-0.37±0.21 


1.20 


-0.04±0.29 


-0.19±0.16 


-0.24±0.14 


-0.30±0.10 


-0.37±0.11 


-0.40±0.07 


-0.35±0.14 


-0.31±0.21 


1.25 


-0.06±0.21 


-0.17±0.12 


-0.20±0.11 


-0.25±0.08 


-0.31±0.09 


-0.34±0.06 


-0.27±0.11 


-0.25±0.18 


1.29 


-0.07±0.14 


-0.14±0.08 


-0.16±0.07 


-0.19±0.05 


-0.23±0.06 


-0.25±0.04 


-0.20±0.08 


-0.18±0.13 


1.34 


-0.06±0.06 


-0.08±0.04 


-0.09±0.03 


-0.11±0.02 


-0.13±0.02 


-0.14±0.02 


-0.10±0.04 


-0.09±0.06 


1.40 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


1.51 


0.25±0.02 


0.24±0.03 


0.24±0.03 


0.24±0.03 


0.25±0.02 


0.26±0.02 


0.18±0.05 


0.16±0.08 


1.66 


0.75±0.15 


0.67±0.08 


0.63±0.08 


0.62±0.06 


0.61±0.07 


0.63±0.07 


0.39±0.12 


0.34±0.15 


1.81 


1.19±0.26 


1.07±0.13 


1.00±0.12 


0.97±0.08 


0.94±0.11 


0.94±0.10 


0.57±0.18 


0.48±0.20 


1.96 


1.42±0.24 


1.37±0.14 


1.27±0.12 


1.23±0.09 


1.20±0.11 


1.18±0.10 


0.71±0.24 


0.55±0.23 


2.11 


1.57±0.23 


1.56±0.17 


1.44±0.15 


1.39±0.12 


1.36±0.15 


1.31±0.10 


0.75±0.33 


0.53±0.30 


2.26 


1.53±0.25 


1.55±0.20 


1.41±0.18 


1.35±0.15 


1.31±0.18 


1.25±0.12 


0.59±0.41 


0.32±0.37 


2.41 


1.21±0.25 


1.26±0.22 


1.10±0.19 


1.04±0.17 


1.00±0.20 


0.92±0.12 


0.22±0.46 


-0.08±0.41 


2.55 


0.75±0.25 


0.80±0.22 


0.64±0.20 


0.58±0.17 


0.54±0.20 


0.46±0.12 


-0.25±0.46 


-0.55±0.41 


2.70 


0.29±0.25 


0.33±0.22 


0.17±0.19 


0.12±0.17 


0.07±0.20 


0.00±0.12 


-0.71±0.46 


-1.01±0.41 


2.85 


-0.14±0.25 


-0.09±0.22 


-0.25±0.19 


-0.31±0.17 


-0.35±0.20 


-0.43±0.12 


-1.14±0.46 


-1.44±0.41 


2.99 


-0.59±0.25 


-0.54±0.22 


-0.70±0.20 


-0.76±0.17 


-0.80±0.20 


-0.88±0.12 


-1.59±0.46 


-1.89±0.41 


3.19 


-1.20±0.25 


-1.16±0.21 


-1.31±0.19 


-1.37±0.17 


-1.41±0.20 


-1.48±0.12 


-2.19±0.45 


-2.47±0.40 


3.53 


-1.83±0.29 


-1.79±0.19 


-1.90±0.18 


-1.94±0.15 


-1.97±0.18 


-2.02±0.13 


-2.68±0.39 


-2.86±0.34 


3.86 


-2.06±0.35 


-2.04±0.18 


-2.07±0.19 


-2.10±0.14 


-2.09±0.17 


-2.12±0.15 


-2.72±0.35 


-2.73±0.37 


4.20 


-2.02±0.41 


-2.01±0.21 


-1.98±0.21 


-1.99±0.15 


-1.96±0.19 


-1.96±0.19 


-2.51±0.38 


-2.37±0.49 


4.54 


-1.83±0.44 


-1.82±0.23 


-1.77±0.22 


-1.77±0.16 


-1.73±0.20 


-1.71±0.21 


-2.24±0.41 


-2.03±0.56 


4.87 


-1.58±0.43 


-1.58±0.23 


-1.53±0.22 


-1.53±0.16 


-1.49±0.20 


-1.48±0.21 


-2.01±0.41 


-1.81±0.55 


5.21 


-1.32±0.43 


-1.31±0.23 


-1.26±0.22 


-1.26±0.16 


-1.23±0.20 


-1.21±0.21 


-1.74±0.41 


-1.54±0.55 



Table 2. Continue. 
























Starburst Galaxies 








log(A) 








log(i 2 5 M m 


1 Lq) bins 








(jim) 


7±1 


8.5 ±0.5 


9.2 ±0.2 


9.6 ±0.2 


10.0 ±0.2 


10.4 ±0.2 


10.8 ±0.2 


11.5± 0.5 


-1.00 


-1.54±0.04 


-2.16±0.51 


-2.45±0.56 


-2.73±0.57 


-2.95±0.46 


-3.20±0.45 


-3.25±0.51 


-3.81±0.61 


-0.76 


-1.26±0.06 


-1.95±0.49 


-2.20±0.53 


-2.48±0.56 


-2.69±0.47 


-2.95±0.46 


-3.01±0.54 


-3.53±0.53 


-0.51 


-1.04±0.04 


-1.63±0.39 


-1.84±0.41 


-2.08±0.42 


-2.29±0.36 


-2.52±0.36 


-2.64±0.42 


-3.12±0.40 


-0.27 


-0.85±0.01 


-1.26±0.30 


-1.44±0.27 


-1.62±0.27 


-1.84±0.26 


-2.03±0.26 


-2.22±0.28 


-2.66±0.32 


-0.03 


-0.65±0.01 


-0.91±0.32 


-1.07±0.21 


-1.19±0.24 


-1.41±0.29 


-1.57±0.25 


-1.83±0.25 


-2.23±0.34 


0.21 


-0.45±0.02 


-0.70±0.36 


-0.82±0.22 


-0.93±0.27 


-1.14±0.31 


-1.28±0.25 


-1.56±0.26 


-1.92±0.33 


0.46 


-0.73±0.02 


-0.93±0.29 


-1.00±0.19 


-1.08±0.23 


-1.25±0.26 


-1.37±0.23 


-1.56±0.23 


-1.93±0.30 


0.64 


-0.94±0.15 


-1.10±0.20 


-1.10±0.21 


-1.17±0.26 


-1.27±0.29 


-1.38±0.31 


-1.45±0.30 


-1.89±0.34 


0.83 


-0.57±0.14 


-0.73±0.19 


-0.74±0.20 


-0.82±0.23 


-0.93±0.26 


-1.03±0.27 


-1.02±0.30 


-1.31±0.37 


0.87 


-0.14±0.15 


-0.32±0.19 


-0.34±0.20 


-0.45±0.23 


-0.60±0.27 


-0.80±0.27 


-0.90±0.25 


-1.31±0.27 


0.91 


-0.25±0.14 


-0.41±0.19 


-0.41±0.20 


-0.47±0.23 


-0.53±0.26 


-0.55±0.27 


-0.53±0.30 


-0.89±0.31 


0.95 


-0.64±0.13 


-0.74±0.16 


-0.71±0.17 


-0.74±0.20 


-0.77±0.22 


-0.76±0.24 


-0.72±0.28 


-1.06±0.30 


0.99 


-0.73±0.14 


-0.87±0.16 


-0.87±0.17 


-0.91±0.19 


-0.94±0.20 


-0.97±0.18 


-1.04±0.17 


-1.31±0.21 


1.04 


-0.50±0.14 


-0.65±0.13 


-0.67±0.15 


-0.73±0.15 


-0.78±0.17 


-0.84±0.15 


-0.92±0.14 


-1.19±0.15 


1.08 


-0.14±0.14 


-0.27±0.14 


-0.27±0.16 


-0.32±0.17 


-0.38±0.19 


-0.40±0.20 


-0.36±0.24 


-0.72±0.25 


1.12 


-0.22±0.13 


-0.32±0.12 


-0.31±0.14 


-0.33±0.15 


-0.35±0.16 


-0.35±0.13 


-0.37±0.12 


-0.56±0.15 


1.16 


-0.35±0.12 


-0.42±0.07 


-0.42±0.10 


-0.44±0.09 


-0.47±0.09 


-0.48±0.05 


-0.49±0.05 


-0.57±0.12 


1.20 


-0.32±0.11 


-0.39±0.07 


-0.39±0.10 


-0.41±0.08 


-0.43±0.08 


-0.41±0.04 


-0.43±0.06 


-0.47±0.11 


1.25 


-0.26±0.08 


-0.32±0.05 


-0.32±0.07 


-0.33±0.06 


-0.36±0.06 


-0.37±0.03 


-0.38±0.03 


-0.41±0.09 


1.29 


-0.21±0.06 


-0.24±0.03 


-0.24±0.05 


-0.25±0.04 


-0.26±0.04 


-0.27±0.02 


-0.28±0.02 


-0.31±0.06 


1.34 


-0.12±0.02 


-0.13±0.01 


-0.13±0.02 


-0.14±0.02 


-0.14±0.02 


-0.15±0.01 


-0.15±0.01 


-0.16±0.03 


1.40 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


O.OOiO.OO 


1.51 


0.29±0.01 


0.27±0.02 


0.28±0.02 


0.28±0.02 


0.28±0.03 


0.29±0.04 


0.29±0.03 


0.31±0.05 


1.66 


0.74±0.03 


0.66±0.07 


0.67±0.07 


0.66±0.08 


0.65±0.08 


0.68±0.10 


0.67±0.08 


0.70±0.11 


1.81 


1.10±0.06 


0.95±0.12 


0.96±0.11 


0.94±0.12 


0.92±0.12 


0.94±0.13 


0.91±0.11 


0.90±0.15 


1.96 


1.28±0.06 


1.09±0.15 


1.10±0.13 


1.07±0.15 


1.04±0.14 


1.04±0.14 


0.99±0.13 


0.90±0.17 


2.11 


1.38±0.05 


1.14±0.19 


1.14±0.17 


1.11±0.19 


1.07±0.17 


1.05±0.17 


0.97±0.17 


0.81±0.22 


2.26 


1.29±0.06 


1.01±0.22 


1.01±0.21 


0.98±0.23 


0.93±0.21 


0.89±0.21 


0.80±0.22 


0.58±0.28 


2.41 


0.95±0.06 


0.65±0.24 


0.65±0.23 


0.61±0.25 


0.56±0.23 


0.51±0.23 


0.41±0.24 


0.15±0.31 


2.55 


0.48±0.06 


0.19±0.24 


0.18±0.24 


0.15±0.25 


0.10±0.23 


0.04±0.23 


-0.06±0.25 


-0.32±0.31 


2.70 


0.02±0.06 


-0.27±0.24 


-0.28±0.23 


-0.32±0.25 


-0.37±0.23 


-0.42±0.23 


-0.52±0.25 


-0.78±0.31 


2.85 


-0.40±0.06 


-0.70±0.24 


-0.70±0.23 


-0.74±0.25 


-0.79±0.23 


-0.85±0.23 


-0.95±0.25 


-1.20±0.31 


2.99 


-0.86±0.06 


-1.15±0.24 


-1.16±0.24 


-1.19±0.25 


-1.24±0.23 


-1.30±0.23 


-1.40±0.25 


-1.66±0.31 


3.19 


-1.47±0.05 


-1.76±0.24 


-1.76±0.23 


-1.79±0.25 


-1.84±0.23 


-1.90±0.23 


-2.00±0.24 


-2.25±0.31 


3.53 


-2.11±0.03 


-2.29±0.22 


-2.28±0.21 


-2.30±0.23 


-2.35±0.21 


-2.40±0.22 


-2.47±0.22 


-2.69±0.31 


3.86 


-2.35±0.00 


-2.39±0.20 


-2.36±0.21 


-2.37±0.23 


-2.42±0.20 


-2.45±0.22 


-2.48±0.21 


-2.65±0.33 


4.20 


-2.32±0.04 


-2.22±0.21 


-2.18±0.23 


-2.18±0.24 


-2.22±0.20 


-2.25±0.24 


-2.24±0.22 


-2.36±0.36 


4.54 


-2.14±0.05 


-1.98±0.22 


-1.93±0.24 


-1.92±0.25 


-1.96±0.21 


-1.99±0.25 


-1.96±0.23 


-2.06±0.38 


4.87 


-1.89±0.05 


-1.74±0.22 


-1.69±0.24 


-1.68±0.25 


-1.73±0.21 


-1.75±0.24 


-1.73±0.23 


-1.83±0.37 


5.21 


-1.62±0.05 


-1.48±0.22 


-1.42±0.24 


-1.42±0.25 


-1.46±0.21 


-1.48±0.24 


-1.46±0.23 


-1.56±0.37 



